Non-Gaussian Signatures from the Postinflationary Early Universe 
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We consider contributions to non-Gaussianity of the Cosmic Microwave Background (CMB) from 
remnants of phase transitions in the very early universe. Such signatures can optimistically be used 
to discover evidence of new particle physics through cosmological observations. More conservatively 
they may provide an obstacle to extracting information about the non-Gaussian nature of primordial 
density fluctuations from any detection in the CMB. We study this explicitly by computing the 
bispectrum from global textures, which occur in a wide class of particle physics models. 
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Measurements of the Cosmic Microwave Background 
(CMB) temperature anisotropy have reached unprece- 
dented precision, and there is now a growing interest [1- 
5] in the non-Gaussianity of the temperature field. Most 
of the current focus is on the exciting possibility of in- 
terpreting any detection as a signature of inflationary, 
or other models of the early universe. To this end con- 
siderable work has been carried out to predict the non- 
Gaussianity expected in the primordial spectrum of per- 
turbations from inflation [6-15] and cyclic [16] models. 

However, as has been discussed previously [17-23], 
primordial effects are not the only way in which non- 
Gaussianity can arise in the observed anisotropy spec- 
trum of the CMB. As photons travel from the surface of 
last scattering, they traverse a variety of structures, many 
of which may imprint secondary non-gaussianities on the 
spectrum. An example of this is provided by the cor- 
relation between the Sunyaev-Zel'dovich effect and weak 
lensing as photons pass through hot galaxy clusters. 

Another possibility is that phase transitions in the 
early universe can introduce a new source of non- 
Gaussianity, which needs to be considered when trying 
to differentiate primordial origins from secondary effects 
in any detected signal. Phase transitions in which topo- 
logical defects form are particularly interesting [24] , since 
these may source non-Gaussianity actively throughout 
cosmic history, allowing for an integrated effect. Thus, 
even though topological defects cannot be the primary 
agent in structure formation, they may contribute a sig- 
nificant non-Gaussianity in the CMB temperature field. 

In this letter we focus on the example of a rather spe- 
cial class of topological defects - textures [25-28]. These 
occur whenever a particle physics theory possesses a 
global symmetry group G that is broken to a subgroup 
H in such a manner that the third homotopy group 7r3 of 
the vacuum manifold Ad = G/H is nontrivial. A simple 
example is given by the Lagrangian density 
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choice of the symmetry breaking potential 
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with A, 7] constants, breaks the symmetry group G — 
50(4) down to H = SO (3), with ir 3 (G/H) ^ Z. 

Textures consist entirely of gradient energy, which 
evolves before ultimately unwinding to topological triv- 
iality in a small region of spacetime, and radiating its 
energy out to infinity. Our goal here is to estimate the 
non-Gaussianity from textures by computing the bispec- 
trum, illustrating the care required when interpreting any 
detected non-Gaussianity as a primordial signal. 

As usual, the fluctuations of the temperature of the 
CMB can be decomposed into spherical harmonics via 
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The bispectrum is then the 3-point correlation function 
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and one observes the angle-averaged bispectrum 
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where <f>i (i = 1,...,4) are real scalar fields, and the 



Here, in order to study the non-Gaussianity that tex- 
tures source, we compute the bispectrum (4), elaborating 
on the analytical results of [28, 29]. 

Any random distribution of scalar field orientations 
will contain a significant amount of gradient energy, 
which will relax to vacuum. Textures themselves are 
entirely gradient energy and, due to topology, a con- 
figuration will collapse until there is sufficient gradient 
energy in a small enough volume that it becomes ener- 
getically favorable for the field to unwind. Both cases 
have an evolving gravitational potential as the config- 
uration reaches triviality by radiating goldstone modes. 
As articulated in [25] , the effect of a particular texture on 
CMB photons depends on whether the photons traverse 
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the texture while it is unwinding or afterwards. Thus, 
textures create hot and cold spots in the CMB. Simula- 
tions [27] , show an analogous effect even if the field con- 
figuration in question does not carry topological charge. 

The defect network, and consequently the network of 
hot and cold spots, is expected to scale, i.e. to look sta- 
tistically the same at any time, with characteristic length 
scale normalized to the horizon size. Following [28], we 
model the spot network with a distribution function 
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where n is the average number of texture-producing con- 
figurations per horizon volume, y(t) = log 2 (to/t) and to 
is the cosmic time today. (Note that we have assumed a 
matter dominated universe, since the most important ef- 
fects come form later textures, during this epoch.) Each 
spot will then contribute a temperature fluctuation 
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where a, and Si are the brightness and the profile, re- 
spectively, of the i-th spot, and 9i is the angle in the sky 
as measured with respect to the center of the spot. A 
spot appearing at time y has a typical angular size 9 s (y) 
which is a fraction d s of the horizon angular size 
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The contribution of the i-th texture spot to a; m is then 
a} m = a i St(y)Y l * m (n i ) , (8) 
where Si(y) is the multipole of the spot profile 



Sl(y) = 2ir J ^(e i ,y)P l (cose i )d(cos0 i ) . (9) 

Summing over the distribution of spots in the sky from 
last scattering until today (and assuming no correlations 
between spots laid down at different times, between spot 
brightness and profile, and that spot brightness doesn't 
depend on the time at which the spot was laid down) one 
finds that the global texture bispectrum takes the form 

B hhh = 7*1/2/3 ( a )^hl 2 l 3 ' ( 10 ) 
where is the integral over the profile multipoles 
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with y is the value of y at last scattering and where 
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FIG. 1: Top panel: equilateral bispectrum for textures, _B ; text 
as a function of I. Bottom panel: S/N 
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Here (• • • ) is the average over simulated spot dis- 
tributions. As we do not have at our disposal the 
full spot brightness distribution from numerical simula- 
tions, we assume that any hot (cold) spot has the same 
Ohot > 0(a co id < 0) corresponding to the averaged maxi- 
mum(minimum) brightness of simulations. Then we may 
write (a 3 ) = x(a 2 } (\a\), and introduce the parameter 
x = (a)/(\a\) (assumed to be ^-independent) encoding 
the asymmetry between cold and hot spots [29]. 

For simplicity we specialize to the gaussian profile, 



Si(y) =2tt[1- cos (9 s (y))} exp 
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and focus on calculating the contribution of textures to 
the equilateral bispectrum (for which li = l 2 = I3 = I, not 
to be confused with the "equilateral shape bispectrum" ) 
on large scales. We choose n = 1, d s — 0.1 (as suggested 
by simulations [27], ensuring that we include the effects 
of all scalar field gradients) and x ~ 0.07 [29]. We leave 
(a 2 ) oc rf as a free parameter related to the symmetry 
breaking scale 77. Note that the observed CMB angular 
power spectrum bounds r\ < 10~ 2 Mp when d s = 0.1, 
(using the C10 normalization from WMAP5 [5]). 

We focus on the angle-averaged bispectrum (normal- 
ized as in [24]) in the equilateral case, h — h = h 
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In fig. 1 we plot the texture equilateral bispectrum B\ cyit , 
finding it to be a monotonically decreasing function of 
multipole index; i.e. textures contribute more non- 
Gaussianity on larger scales (in agreement with [28], in 
which the angular spectrum Ci was studied.). In the 
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FIG. 2: The total S/N for the equilateral bispectrum Bffi* of 
textures as a function of the symmetry breaking scale rj. 

same figure we plot the theoretical signal-to-noise ratio 
S/N for the equilateral bispectrum Bm, which is also a 
monotonically decreasing function of I. For the variance 
of Bui, we use the equilateral version of an expression 
in [30], valid for mildly non-Gaussian distributions, 

a? = (B? u )-(B m ) 2 *6Ci , (15) 

where C; includes the power spectrum of the detector 
noise, calculated using the analytical expression of [31], 
and we have included the contribution from textures. 
Therefore we have (S/N) [ w \B m /Cf /2 \. In fig. 2 we 
plot the total S/N as a function of n. As the symmetry 
breaking scale increases, the signal increases at a bigger 
rate than the associated noise, therefore the ratio S/N is 
found to be an increasing function of 77, with S/N ~ 0(1) 
for 77 <~ 5 • 10~ 3 Mp. To gain some physical insight, wc 
may look at the angle-averaged bispectrum density for the 
equilateral case: 

Br\y) = V^7ih I (« 3 > N(y) Si(yf . (16) 

It is worth noticing that the geometrical factor (12) en- 
ters eq. (16) quadratically, therefore there are no oscilla- 
tions associated with the alternating sign of the 3j-wigncr 
symbol. In the case of gaussian profiles this reads 

Br\y) = 32^/^<a 2 > 3/2 ^m(2) ^ - if 

•[l-cos(0 s (y))] 3 e -H^)X 2 fe) t (17) 

and in fig. 3 we plot B\ ext (y) as a function of the "time" 
variable y. The width of the curves increases as a func- 
tion of the multipolc index I, while the height of the peaks 
decreases as a function of I. The tilt of the bispectrum 
in fig. 1 depends on the balance of these two effects. It is 
worth pointing out that, as I decreases, not only do the 
curves narrow, but the peaks shift to smaller values of 
the time y; i.e. lower multipoles receive most contribu- 
tion from the most recently formed textures (A similar 
effect was found in [28] for the power spectrum itself.) 
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FIG. 3: The texture equilateral bispectrum density, _B* ext (j/), 
for d s = 0.1, n = 2 • 1CT 3 Mp, and from top to bottom 4 < 
I < 16. The curves have been normalized to b4,(y max ) 
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FIG. 4: Equilateral bispectrum density for I = 4, B\ ext (y), 
for three different values of the density parameter; d s =0.1 
(solid line), d s — 0.2 (dashed line) and d s = 0.4 (dotted line) 
(normalized to the case d s = 0.1). In all cases n = 2-10~ 3 M P . 



It is interesting to study B\ cxi function of the 

impact parameter d s at fixed multipole I. We do this 
in fig. 4 for I = 4 for three different values of the impact 
parameter. As can be seen, the position of the peak shifts 
to higher values of y for increasing values of d s . Thus, the 
larger the value of d s , the earlier the time at which the 
main contribution to the non-Gaussianity is generated. 
It is worth noting that, as d s is changed, the position 
of the peak changes such as to keep the angular size of 
the spot nearly fixed (for fixed I). This is as one might 
expect, since the spots are uncorrelated. Thus, the main 
contribution to non-Gaussianity at a given multipolc I 
arises at times at which the textures have an angular size 
of order the angular size of the multipolc ~ 4tt /l(l + 1). 

Post-inflationary early universe physics may generate 
non-Gaussian signatures in the CMB. Such signatures 
may be useful in identifying new physics in upcoming 
missions, but it is also important to understand them in 
order to be able to interpret any signal as evidence for 
primordial physics. As a concrete example we have com- 
puted the cosmic texture bispectrum, in a way differing 
from previous analyses [28, 29, 32-38] in the analytical 
modeling and the use of updated datasets and current 
bounds on the allowed texture energy scale. Our result 
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FIG. 5: Ratio of the bispectrum from textures, for r\ — 
7 ■ 1CT 3 M P , to the bispectrum from DBI inflation [39], with 
speed parameter 7 = 30, as a function of the multipole index. 



may be compared to the primordial bispectrum from, 
say, DBI inflation, which is of the equilateral form [39]. 
In fig. 5 we show an example of the ratio of these bispec- 
tra, which is a nontrivial function of the multipole index 
I, since the individual bispectra have different shapes. 
Observational constraints [5] are commonly quoted in 
terms of the quantity /nl- We may estimate the equi- 
lateral /nl from textures by comparing Bf cxt to that ob- 
tained from the 3-point correlation function of the cur- 
vature perturbation 

= (-jj/NLAf) A -^0 

1 (18) 

evaluated on an equilateral triangle k± — k 2 — k 3 . The 
resulting estimate is a function of the symmetry breaking 
scale, with |/^ text | ^ 5 • 10 14 (r//M P ) 6 with a value of 
O(10) for r\ ~ 5 • 10" 3 M P . 
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